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Abstract 

We have investigated the effects of primeval starbursts on the galactic structure, showing that the dy- 
l/^ namical structure of young galaxies is changed by relaxation of a stellar system formed by a starburst-driven 
^\ outflowing gas (superwind). When a superwind interacts with halo gas, an expanding dense supershell of 
shocked gas is formed. Using a similarity solution for the expanding shell and a condition for the gravita- 
tional instability of a gaseous shell, we show that an expanding gaseous shell with a mass of several 10"'^ "^Mq 
^ and a radius of several kpc becomes gravitationally unstable to form stars. A stellar shell is thus formed, 
relaxes and will evolve into a fat stellar system. 

In order to investigate the fate of the stellar shell and its dynamical influence on the host galaxy, we 
performed three dimensional #-body simulations of a stellar shell plus disk system. The evolution of the 
stellar shell was computed under two types of environment: an external flattened potential and a live stellar 
disk. We find that the relaxation process and resultant structure are significantly affected by the external 
disk potential, and also depend on the velocity dispersion of the initial shell stars and the mass of the shell. 
The final galactic structure ranges from a thick disk and a compact bulge with a high density core to a thin 
disk with a diffuse bulge. Isodensity maps of the computed disk plus bulge system resemble isophotes of SO 
galaxies. We suggest that the effects of primeval starbursts would explain some observational properties of 
SO galaxies. 
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Ol. Introduction ing part. This model has long been the basis for under- 

standing bulge formation, although there have appeared 

Observations have suggested that a substantial frac- several ideas that concern bulge formation after galaxy 
tion of galaxies underwent starbursts (Soifer et al.l986). ^isks formed. Combes et al. (1990) showed that a bar- 
It is plausible that when the universe was younger, the potential m the central region excites vertical motions 
number density of galaxies was high and that galaxy in- ^^ars to form a boxy bulge; Pfenniger and Norman 
teraction was frequent. This would have triggered rapid (1990) showed that contraction of disk stars toward the 
gas fueling into galaxy nuclei, and starbursts would hkely center m a barred potential resuKs m an increase m ve- 
have occurred in a young and gas-rich disk galaxy. It is, lo^ity dispersion perpendicular to the disk, and yields 
therefore, worthwhile to investigate how primeval star- ^ bulge-like distribution. However, m these studies ef- 
bursts afl^ected the evolution of structure and properties ^^cts of starbursts induced by galaxy-galaxy interactions 
of young galaxies young galaxies were not considered. Recently, Sofue 

The formation of a galactic bulge in a disk galaxy is ^nd Habe (1992) suggested formation of a central bulge 

one of the long-standing questions in the evolution of ^y a primeval starburst, and showed that stars formed 

galaxies. Larson (1976) showed that a bulge formed dur- starburst-driven gases would evolve into a bulge sys- 

ing an initial spherical contraction of a central core in ^^m, whose size depends on the strengths of starburst 

a protogalactic gas sphere. This was foUowed by disk activity. However, their model is rather speculative and 

formation from a pancake contraction of the outer rotat- quahtative. They did not consider self-gravity of a newly 
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formed stellar system and gravitational interaction of the 
bulge and disk system. In the bulge formation hypothe- 
sis by starbursts, these effects should be considered more 
carefully. 

In the present paper, we intend to show that a substan- 
tial fraction of bulge stars could have originated in star 
formation in a starburst-driven supershell in primeval 
galaxies. In section 2, we evaluate properties of a stellar 
shell system by using a similarity solution of an expand- 
ing gas shell in an accreting ambient gas and a condition 
for gravitational instability. To investigate the evolution 
of this stellar system under the influence of the non- 
spherical external potential of a galactic disk, we perform 
three dimensional N-body simulations. Our methods and 
models are described in section 3. The models are cho- 
sen to show the effects of the initial conditions on the 
final structure and on the relaxation process. The initial 
conditions are within the limits provided by our analyt- 
ical calculations. Numerical results are shown in section 
4. Finally in section 5, we summarize our results and 
discuss implications of our results on the origin of SO 
galaxies. 

2. Gaseous and Stellar Shell Formation by a Su- 
perwind 

2.1. Superwmd Shell 

Many hydrodynamical simulations have shown that 
a superwind produced by multiple supernovae interacts 
with the ambient gas and forms a hot cavity and a cooled 
dense gas shell (e.g. Tenorio-Tagle, Bodenheimer, Rozy- 
cka 1987; Tomisaka, Ikeuchi 1988; Mac Low, McCray 
1989). Star formation in the gaseous shell would occur, if 
the shell is radiatively cooled and becomes gravitationally 
unstable (McCray, Kafatos 1987, Ostriker, Cowie 1981). 

Recently, some high redshift radio-galaxies were found 
to have analogous properties to nearby far infrared galax- 
ies (McCarthy, Persson, West 1992), in and around which 
intense star formation is ongoing. In these candidates of 
primeval galaxies, a superwind and a secondarily-formed 
stellar system along the radio-lobe axis have been ob- 
served (MacCarthy et al. 1987). It is suggested that an 
expanding super bubble in a gaseous halo triggers forma- 
tion of a stellar system (Heckman, Armus, Miley 1990). 
Therefore, it seems reasonable to suppose that the stel- 
lar systems observed around high-redshift galaxies were 
formed by these processes. 

We first estimate properties of a gaseous shell that 
would be formed by an interaction of a superwind with 
an accreting ambient gas. Since primeval starbursts and 
their interaction with proto-galactic gas sphere have not 
been fully studied, we here assume that they have similar 
properties to those observed in nearby active starburst 
galaxies. 



From observations of double emission-lines with line 
splitting of 200 — 600 km/s in nearby far infrared galaxies, 
the presence of a "superwind" has been suggested with 
its "wind luminosity" of the order of _Lw > IQ^'^Lq ~ 
4 X 10^5 erg s"! (e.g. NGC 3690, IC 649, Arp 220, and 
NGC 6240), to ~ 10^^1.0 in IRAS 00182-^7112 (Heck- 
man, Armus & Miley 1990). Interaction of such a su- 
perwind with the ambient intergalactic/halo gas would 
produce a dense shell of shocked gas (superwind shell). 
It has been suggested that a stellar shell would be formed 
from such a dense gas shell, and that this could ex- 
plain the origin of stellar shells around E/SO galaxies 
(Umemura, Ikeuchi 1987; Williams, Christiansen 1985; 
Fabian, Nulsen, Stewart 1980), although mergers are the 
commonly accepted explanation of the shells and the rip- 
ples (e.g. Prieur 1990). A smaller scale but more massive 
stellar shell could form in the halo nearer to the nucleus. 

Umemura and Ikeuchi (1987) presented a similarity so- 
lution of an expanding gaseous shell produced by a galac- 
tic wind interacting with an accreting ambient gas in a 
spherical galactic gravitational field. Their similarity so- 
lution of the expanding shell is given by rs{t) = (at)^^^, 
when the initial density distribution of the accreting am- 
bient gas is Pg{r) = Pc{r/rc)~^^^ , where pc and Tc are a 
central density and a core radius, respectively. The crit- 
ical wind luminosity -Lm,crit , above which the wind can 
blow outward against the ram pressure of accreting gas, 
is given by 

X — --s — erg s \ (1) 

when a = 0.21(kpc)~^''^(10® yr)""'^, where ric is the num- 
ber density of the gas at the galactic center and Mg is a 
mass of the galaxy. 

If the energy released from the starburst region were 
less than -Lw,crit) a super bubble confined in a dense 
gas shell with zero mean velocity would be formed (e.g. 
Tenorio-Tagle, Bodenheimer, Rozycka 1987; Tomisaka, 
Ikeuchi 1988; Mac Low, McCray 1989). It is most in- 
teresting when the ejected energy, from a starburst 
region is comparable to -Lw,crit) because in this case the 
shell attains maximum radius and a large amount of am- 
bient gas can be swept up into it. 

In this case, we can estimate the shell radius fg, and 
mass Ms, by using the similarity solution with the above 
condition for a, as 

r-s(i) = 7.8 (^) ' kpc, (2) 
and 

MM) = l.2xW'°( ^^if^^^V 

^ ' Vo.01cm-3y VlOkpc J 
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Mr 



(3) 



2.2. Formation of a Stellar Shell 

Next, we estimate a time when the expanding gas shell 
begins to fragment due to its self-gravity. Kinetic, grav- 
itational, and thermal energy of a small part of the thin 
gas shell between an angle 6 and 6+A6, the mass of which 
is S(rs)(rsA0)^ , with S(rs) being the surface density of 
the shell at radius fg, are defined as 



Ei,(t) 
and 



1 



S(r,)(r,Aef(r,Aey 

GS{r,f{r,Aey 
" r,A0 



S(r,)(r,Aef 



IJ-muij - 1) 



(4) 
(5) 

(6) 



respectively, where a is a constant of order unity, mpj 
is the mass of a hydrogen atom, /j, = 0.61 is the mean 
molecular weight for a fully-ionized gas of #He/^H = 0.1, 
')/ = 5/3 is the adiabatic exponent, and Tg is a temper- 
ature of the gas shell. Using the adiabatic strong shock 
condition, the gas temperature Ts(t) of the shell is 



Ts(t) 



3iJ,miiVs(ty 
I6k 



(7) 



where Vs(t) = rs(i)- The total energy of the gas shell 
E^otai = E]^ + Eg + becomes less than zero, when 



(aGSrs)^ - 2rl 



IJ-muij - 1) 



> 0. 



(8) 



Hence, equation (8) is a condition that the gas shell be- 
comes gravitationally unstable. 

Let us turn to a cooling time of the shocked gas at the 
leading edge of the expanding gas shell. From the energy 
conservation law for radiatively cooling monatomic gas 
with initial temperature Tg, the cooling time tcooi of the 
shocked gas is given by 



^cool 



lOfc^T. 



2-S 



3(2 - 6)Apg(r,)v, 



2.7 X lO-' 



Tlr 



108 yr 



0.01 cm-3 

0.73 



10 kpc 



(9) 



where we assume the cooling function A = AT^ , with 
6 = -0.6, and A = 6.2 x lO-^^erg cm^ s"! (McKee, 
Cowie 1977). Equation (9) shows that tcooi is al- 
ways less than the age of the shell, after t = 2.3 x 
10^ (nc/0.01cm-3)-0 58 (r JlO kpc)-° ^9 yr. Therefore, 



the temperature of the shell quickly cools to ~ lO'* K at 
which it remains because of a secondary heating source: 
photoionization due to the nonthermal continuum emis- 
sion which accompanies the wind as long as the nucleus 
remains active (Williams, Christiansen 1985). 

The critical time, tcrit, when the gas shell just becomes 
gravitationally unstable, is derived from equations (2) 
and (8) with = lO^* K and a = 1, 



tr 



1.6 X lO"" 



0.01 cm-3 



10 kpc 



(10) 



This time scale is comparable to a typical starburst life 
time. The mass and the radius of the shell at tcrit are 



Msitcrh) 



2.6 X 10' 



0.01 cm-3 



10 kpc 



Me, (11) 



(^crit ) 



10.7 



0.01 cm-3 



10 kpc 



kpc, (12) 



respectively. Soon after the starburst is over, the tem- 
perature of the shocked gas immediately cools to several 
10^ K. At this stage the thickness of the gas shell is much 
less than its radius, because the internal pressure of the 
shell cannot support it against the external ram pressure. 

Note that in above discussion we assumed that the 
cooling is very effective. However, if some heating mech- 
anisms exist in the process, the temperature Tg would 
remain higher than lO'* K and the critical time would be 
longer than that given by equation (10). Therefore in a 
higher density halo, a more massive gas shell than that 
estimated here would be expected. 

At t ~ ^crit, fragmentation of the shell occurs due to 
gravitational instability, and star formation is triggered. 
In order to estimate mass of a stellar shell formed from 
the gas shell, we need the star formation efficiency e^. 
It has been suggested that 0.3 < < 0.5 is needed to 
form bound star clusters from molecular clouds (Lada, 
Margulis, Derborn 1985). If we assume in our primeval 
starburst shell was comparable to this Galactic value, 
a thin stellar shell of Mg star ~ several xIO-'^'^Mq with 
f'sjstar ~ several kpc can be formed. 

Further evolution of this stellar shell depends on the 
initial mass function of stars (IMF). If we assume that 
the IMF was similar to that obtained for the solar neigh- 
borhood (Miller, Scalo 1979), the birth rate of massive 
stars which would explode as SN can be estimated using 
the IMF: ^(logM) cx Af-^ ^^ for 10 > M/Mq > 8, and 
for M/Mq > 10. This gives a rate of SN ex- 
plosions as 4 ~ 10 X 10"-^^ yr~-^ pc~^, or 0.05 ~ 0.1 
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yr ^ in a shell of a radius 10 kpc. This rate would 
be maintained for the Jeans time of the gaseous shell, 
Tff ~ 737r/32Gpsheii ~ 10^ yr (Ar, = 100 pc, = 
2.6 X 10"'^'^ Mq, r-g = 10 kpc). Hence, the number of 
SN which occurred in the shell is estimated to be 5 ~ 
10 X 10^ and the total energy release is about 10 erg, 
which is comparable to the total self-gravitational energy 
of the gas sheU, GM^/2rs ~ 3 x 10®'' erg. This energy 
is not enough to destroy the gaseous shell, since about a 
few % of the energy released by SNs changes to kinetic 
energy of SN remnants (Spitzer 1978). Furthermore, if 
the fragmentation of proto-stellar cloud is hierarchical, 
there would be a sharper cutoff for the massive star for- 
mation. Larson (1982) has pointed out that massive stars 
are formed only in high-mass clouds. The Jeans mass in 
the gas sheU is Mj = {irkT / iiGfl"^ p~^l^^ ~ 4 x WMq 
for T = 100 K. According to Ashman (1990), the largest 
fragment in a proto-stellar cloud, whose density fluctua- 
tion is ^(m) cx m.^~^ , is given by 

( B \^ (m^\^ 

x(|^)"*M,. (13) 

Using Mc = Mj, B = -2, and taking nimm = 0.004 Mq 
(Palla, Salpeter and Stabler 1983), we obtain nimax = 
9Mq. For B = —3, we obtain nimax = 0.75 Mq. Silk 
(1977) suggested B = -2.9 for M > 0.6 Mq, if dynami- 
cal dissipation caused by turbulent motion in the proto- 
stellar cloud is effective, which would be likely excited by 
a proto-stellar wind as well as by a hydrodynamical insta- 
bility in the shell. Therefore, it is probable that the IMF 
in the gaseous shell would have a sharp cutoff at around 
~ 1 Mq, and low mass stars dominate in the shell. 

3. Models and Numerical Method 

Suppose that a massive stellar shell comprising small- 
mass stars is formed. Then, our next interest is to know 
how the stellar shell and the disk system evolve. The 
stellar shell would finally relax and evolve into a spherical 
system distributed around the nucleus. We may thus 
expect that the system would evolve into a bulge-like 
structure. In order to investigate how the stellar shell 
would relax, and what shape and kinematic of the final 
"bulge" would exhibit, we perform three-dimensional N- 
body simulations. 

3.1. Model of stellar shells 

The parameters for the stellar shell can only be roughly 
estimated by the preceding analysis. We therefore choose 
a variety of allowed initial conditions, and concentrate 
ourselves on how they affect the final fate of the shell. 



Table 1.. Parameters of the fixed potential models 



RESULT 

Model Ms/Md Figure Type 



Ai 


0.5 


10.0 






III 


A2 


0.5 


5.0 


Fig. 
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III 


A3 


0.5 


1.0 






II 


A4 


0.5 


0.0 






II 


B 


0.1 


1.0 


Fig. 


2 


II 


C 


0.01 


1.0 


Fig. 


1 


I 


D 


0.0 


1.0 






I 



Ms/M<j: the ratio of the mass of the stellar shell and 
the disk mass. ag'. initial velocity dispersion of the shell 
stars. Model D is a test particle model. Type I, II, and 
III mean relaxation types defined in Section 4.1. See 
also Fig. 4. 



We assume that the initial stellar shell is spherically sym- 
metric and thin (Arg/rg = 0.05), and has zero mean ve- 
locity (i.e. no rotation and no radial motion). To model 
a stellar shell, 5000 collisionless particles are randomly 
distributed on a spherical shell. The main parameters of 
the models are the mass of the shell, Mg, and the initial 
velocity dispersion of the shell stars, (Tg. We simulate 
the evolution of the stellar shells in two types of model 
disk: (a) a fixed (i.e. time-independent) axisymmetric 
fiat potential, and (b) a live thin stellar disk modeled by 
#-body selfgravitating particles. 

3.2. Fixed potential models 

The fixed disk potential is assumed to be the 
Miyamoto-Nagai potential (Miyamoto & Nagai 1975): 

^mn(R, z) = = , (14) 

R-^ + (a + vG^Th^y 

where a = 20 kpc, h = 2 kpc, and Mj = l.Ox 10"Mq. In 
this paper we use the following units: [L] = kpc, [M] = 
Mq, [T] = lO^yr, and [V] = [L/T] ~ 10 km s"! . The 
initial radius of the shells is 5.0. 

In Table 1 we summarize the parameters of the initial 
shells in the fixed potential models. Model D is a test- 
particle model in which we do not take into account the 
self-gravity of the shell. 

Since the model A series has a large Mg/M<j, the as- 
sumption of a fixed disk potential may be fiawed. How- 
ever, as we describe in section 4, the evolution of these 
massive shells is similar in both the fixed and lived po- 
tential. 
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Table 2.. Parameters of the live disk models 



Model 




Ns 




OS 


Fi^ 


jure 


SI 


10000 


2000 


0.2 


1.0 






S2 


10000 


5000 


0.5 


1.0 


Fi{ 


5 


S3 


10000 


5000 


0.5 


10.0 
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#<j and Ns are the number of the disk stars and of shell 
stars, respectively. 

Shell of model S2 and that of S3 correspond to A3 and 
Ai of the fixed potential model, respectively. 



3.3. Live disk models 

The second group of models is calculated in order to in- 
vestigate gravitational interaction between the shell and 
disk stars. Both the shell and disk are represented by col- 
lisionless particles (Table 2). The initial density profile 
of the disk is 



Pd(R,z) 



2Trh Ro 



z 



(15) 



where So is a surface density at i? = 0, and Rq and h are 
the radial and disk scale lengths, respectively. 

The vertical scale height h and the vertical velocity dis- 
persion (Tz are determined so that the disk is quasi static 
in the vertical direction. The radial velocity dispersion at 
R = Rq is determined from (TR(i?o) = 3.36(5GSo/k(-Ro), 
where Q is Toomre's disk stability parameter (Toomre 
1963), and we take Q = 1.2. The vertical velocity disper- 
sion (Tz is assumed to be = 0.5(Tr. The scale height 
of the disk is assumed to be constant in the disk, and is 
derived from the hydrostatic condition and the Poisson's 
equation of a very fiattened system. 



1 dpd(Tz 

dz 



Pd 
and 

dz'^ 



5$ 

dz ' 



1 dRF^ 
R dR 



(16) 



(17) 



where _Fr is a radial force of gravity of the disk (Binney, 
Tremaine 1987). 

We also adopt a fixed dark halo potential: ^dark = 
— GMdh/vf^+JV^) where M<jh and ro are the mass and 
core radius of the dark halo. We take ro = 10 and M<jh = 
5.0 M<j = 5.0 X 10"'^"'^ Mq so as to satisfy the criterion of 
a stable disk against the bar instability (Ostriker and 
Peebles 1973). The initial radius of the shells in these 
models is 3.0. 

After several rotations of the disk, a shell is added to 
the disk, and we calculate the time evolution of the shell 
plus disk system. 



3.1 
TV- 
wit h 
are e 



and 



Numerical Method 
■body calculations were performed on a workstation 
GRAPE-3. The basic equations of the calculations 
;iven by 

(18) 



minij (xj — Xi ) 

(^..2 + ,2)3/2 ' 



(19) 



where Xi and m,- are the position and the mass of the i-th 
particle, is the force exerted on it, rij is the distance 
between the i-th and j-th particles, G is the gravitational 
constant, and e is a softening parameter. We adopt a 
constant e = 0.08 and all particles have the same mass 
nii = m in our simulations. Instead of adopting methods 
to reduce the operations such as the particle-mesh 
method with the Fast Fourier Transform or tree algo- 
rithms, the force calculation was performed directly on 
GRAPE. 

GRAPE (GRAvity PipE) is a series of special-purpose 
computers to accelerate solving the above equations, de- 
veloped at the Department of Earth Science and As- 
tronomy, University of Tokyo (Sugimoto et al. 1990). 
GRAPE-3 is a parallel-pipeline system, which attains 
15 Gfiops at peak speed (Okumura at al. 1991). Er- 
ror analysis of the GRAPE was done by Makino, Ito 
and Ebisuzaki (1990): Integration error of the total en- 
ergy per crossing time is about 0.01 percent for a stable 
Plummer sphere using 10^ particles. In our calculations 
the energy is conserved to within 1 percent per dynam- 
ical time. We also performed a number of comparison 
simulations by another numerical method (Particle-Mesh 
method with EFT), and found the same results as those 
obtained using GRAPE. 

All other calculations, such as the time integration, 
were performed on the host computer (workstation). The 
time integration was performed using a standard leap- 
frog method. 

4. Numerical Results 

4.1. Fixed potential models 

The violent relaxation for collisionless N-body systems 
with no disk potential is governed by the initial ratio 
of kinetic to potential energy (e.g. van Albada 1982). 
Relaxation processes of 'cold models' (where this ratio is 
low) is more violent than that of 'hot models'. By adding 
a disk potential, the variety of relaxation processes is 
increased, with (Tg a key determinant. We found that the 
relaxation process of seven models shown in Table 1 can 
be classified into three types; Type I (model C and D), 
Type II (A3 , A4 , and B), and Type III (Ai and A2 ). 



6 



K. Wada 



[Vol. 47, 



The final structure and amount of interaction with the 
disk vary by type. Type III models have larger (Tg than 
Type I or Type II . Type I models are dynamically hotter 
than Type II . 

Figure 1 shows the time evolution of the positions of 
#-body particles projected onto the x-z plane (edge-on 
view) for model C. Due to the disk potential, stars near 
the disk fall toward the disk plane. A dense ring begins to 
form (t = 0.68) and expands in the z- direction (t = 1.28). 
Att = 1.85 a similar ring begins to form in a cigar-shaped 
shell at the disk plane. These processes are repeated and, 
as a result, a multiple-cylinder structure is formed (t = 
3.67 ~ 5.37). The shell stars are finally relaxed to form 
a vague bulge (t = 11.2). The same structure is observed 
even in the non-selfgravitating model D. The multiple- 
cylinder structure is, therefore, a pattern caused by phase 
differences between particles oscillated in a non-spherical 
external potential; phase-mixing relaxes the system in 
model C, as in model D . The self-gravity of the shell is 
not important for the relaxation process in Type I model. 
The time needed for the shell stars to be relaxed in Type 

I models is more than 10. 

Next, we show a typical example of a Type II model, 
in which shells are more massive and dynamically colder 
than in Type I models (Fig. 2; model B). The relax- 
ation process and the final structure are very different 
from those of Type I . In the early stage, gravitational 
instability occurs in the collapsing shell. After the col- 
lapse, the shell expands non-isotropically, and forms a 
torus (t = 0.77 ~ 1.15). Gravitational instability grows 
in the torus, and two clumps are formed in the disk plane 
(t = 1.67). The clumps finally merge at t ~ 3.5 and form 
a compact core with a diffuse halo (t = 8.76). The effect 
of the disk potential is also important for the relaxation 
process in Type II models, although the self-gravity dom- 
inates the process. 

Evolution of models with velocity dispersion larger 
than that in Type II is classified as Type III . The 
cylindrical structure no longer appears in these models 
(Fig. 3). This shows that the disk potential does not af- 
fect the relaxation process so much. In model Ai , a shell 
is relaxed within ~ 2 x 10® yr, and a diffuse prolate bulge 
is formed along the disk plane. 

The square of the initial velocity dispersion (Tg of each 
model is plotted against Mg/M^ in Fig. 4. Three evolu- 
tional type in our results can be explained, if there are 
a critical value of (Tg and Ms/M^- For large Mg/M^, the 
shell evolves as if there is no disk potential, provided 
that (Tg is small enough. This is categorized as Type 

II , and selfgravity of the shell dominates the evolution. 
When Mg/Md and (Tg is small, the shell is strongly af- 
fected by the disk potential and Mg/M^ ~ 0.03 is a crit- 
ical ratio. In a non-spherical potential, a phase-mixing 
occurs during the shell relaxation near the disk plane. 
The vertical velocity dispersion (ctj) of the disk stars in 



the Miyamoto-Nagai potential (eq. (14) to (16)), is crit- 
ical for the phase-mixing by the disk potential to occur. 
This velocity dispersion is shown in Fig. 4 for two radii, 
R = 1.0 and R = 10.0. If the initial velocity dispersion 
of the shell stars is far greater than (T<j, the disk potential 
hardly changes the random motion of shell stars, that is, 
relaxation of the shell stars is not infiuenced by the disk 
potential. On the other hand, if the velocity dispersion is 
smaller than (T<j, the disk potential causes phase-mixing 
of shell stars. Thus evolution process in a disk poten- 
tial depends on the initial velocity dispersion of the shell 
stars ((Tg) as well as the mass of the shell (Mg), therefore 
three types of evolution processes appear in our models. 

The final structure of the remnant of the shell stars 
(bulge) is different among types I, II, or III. If the mass 
of the shell is very small, as in Type I , the bulge has a 
very vague shape (Fig. 1). On the other hand, the bulge 
in Type II (Fig. 2) or in Type III (Fig. 3) has a core- 
halo structure, although the core of Type II is smaller 
than that of Type III . The density profile of the bulge 
is p(r) oc at small r and p(r) oc in the outer 
region. This profile is characteristic in remnant formed 
by the violent relaxation of a coUisionless self-gravitating 
#-body system (e.g. van Albada 1982). 

4.2. Live disk models 

Concerning the shell evolution, results of models SI, 
S2 and S3 are similar to those of the fixed potential case. 

We can see the following points from Fig. 5, which 
shows evolution of (a) disk stars, (b) shell stars and (c) 
both disk and shell stars, in model S2, in which Mg/M^ 
and (Tg are equal to those in model A3 . Shell stars evolve 
in almost the same way as in Type II of the fixed poten- 
tial models; After central concentration, the system of 
shell stars expands and changes into a cylindrical shape. 
In this torus two clumps grow early (t = 0.21) near to 
the disk plane. These two clumps finally merge and settle 
into a dense core with a diffuse halo at t ~ 0.5. 

Disk stars are strongly infiuenced by the shell relax- 
ation; stars in the central region follow a similar evolu- 
tion to the shell stars: a central dense core is formed at 
the central concentration stage of the shell stars, and a 
clumpy structure is also formed at the same position as 
the shell stars (t = 0.21 and 0.31). The disk stars are 
heated up in these processes, and the stellar disk rapidly 
thickens during the relaxation of the shell (see below). 
Moreover it is found that an oval ring density wave ex- 
pands in the disk stars(t = 0.21 ~ 0.5). As a result, the 
spiral arms seen in the initial disk disappear quickly. 

Evolution of model SI, in which the mass of the shell 
is smaller than that of S2 , is similar to the results of 
model S2, although the density wave formed in the disk 
is weak. Model S3, in which the shell has same Mg/M^ 
as in model Ai , evolves differently from model S2 (Fig. 
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6). The density wave generated in the disk is not oval 
and the shell stars do not form clumps as seen in Type 
II models and in S2 model. 

Figure 7 shows the time variation of an average scale 
height of the stellar disk in models SI, S2, S3, and NS. 
Model NS is for comparison, in which a spherical fixed 
potential $ = GMNs/\/f^ + r'^ is added to the live disk 
instead of a stellar shell. The core radius = 0.3 and 
Mns is determined so that the potential energy within 
equals that of the initial shell of model S2. Figure 7 shows 
that a thick disk is formed in models with shells which 
are massive and have a small initial velocity dispersion. 

Figure 8 shows the time evolution of mean vertical ve- 
locity dispersion of the whole disk in model S2, and S3. 
As with the scale height (Fig. 7), the velocity dispersion 
increases rapidly for S2 in an early time. In model S2 
the vertical velocity dispersion in the disk increases most 
rapidly and the disk stars are heated up, where the dense 
clumps are formed (cf. Fig. 5). Since these clumps finally 
merge at the central region of the disk, the hot regions 
also merge, and a region with higher velocity dispersion 
is formed. This results in an increase in the vertical ve- 
locity dispersion and therefore in the scale height in the 
inner region of the disk. The increase of velocity disper- 
sion is caused by a non-linear growth of the perturbation 
in disk stars driven by the evolution of shell stars. Due 
to this non-linear growth, potential energy of the region 
increases, and this causes the heating. Since hot shells 
do not produce clumps in a disk plane, it is reasonable 
that heating up of disk stars is more effective for a dy- 
namically cold shell, i.e. S2 models (Type II ) than for 
S3 (Type III ). 

In Fig. 9, we plot time variation of the angular mo- 
mentum of the shell and disk: About 1.2 % (SI), 5.8 
% (S2), and 3.0 % (S3) of the total disk angular mo- 
mentum is transferred from the disk to the shell during 
relaxation. Due to this angular momentum transfer, the 
bulge formed from the non-rotating shell begins to rotate 
in the same direction of the disk. The mean rotational 
velocity of the bulge is about 5 at i? = 1 in model S2. The 
angular momentum transfer is caused by dynamical fric- 
tion between the clumps and disk stars. When the shell 
is dynamically cold, such as that in model S2, a greater 
number of shell stars evolve into clumps in the disk plane 
(Fig. 5). In fact, the angular momentum transfer is found 
to be highest in model S2. 

The bulges of model S2 and S3 correspond to that of 
Type II and Type III , respectively. Figure 10 shows final 
radial density profiles of the shell stars for model S2 and 
S3. The bulge of both models have a core {p{r) cx r~^) 
and halo {p{r) cx r~'*) structure. However the radius of 
the core in the bulge of model S2 is one third of that in 
model S3, and the central density of model S2 is larger 
than that of model S3. Makino, Akiyama & Sugimoto 
(1990) showed that the surface brightness distribution of 
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a stellar system with density distribution proportional to 
r~'* in the outer region and in the inner region can be 
fit by any de Vaucouleurs type r^l"^ law (de Vaucouleurs 
1948). 

5. Discussion 

5.1. Summary 

We have shown that starbursts in galaxies with an ac- 
creting gas halo can produce a stellar shell with a mass of 
several IQ^''-' Mq and a size of several kpc by using a simi- 
larity solution and a condition for gravitational instabil- 
ity of a gaseous shell (Section 2). In order to investigate 
the evolution of the shell and its dynamical effects on 
the structure of the host galaxy, we performed three di- 
mensional #-body simulations. The initial stellar shells 
were simply assumed to be thin and spherical, and were 
evolved in two case of disk models: a fixed disk poten- 
tial and a live stellar disk. We used the latter model to 
calculate the structural evolution of the disk caused by 
the gravitational interaction of the shell and disk stars. 
The main parameters are the initial velocity dispersion of 
the shell stars and the mass of the shell. Our numerical 
results are summarized as follows. 

1. The relaxation process and resultant structure are 
significantly affected by the external potential, and 
also depend on the velocity dispersion of the initial 
shell stars as well as the mass of the shell. We found 
two new types of relaxation which have not been 
known in the previous collapse simulation of colli- 
sionless system without external disk potential; (a) 
Relaxation is mainly due to phase-mixing caused by 
the disk potential. During the relaxation multiple- 
cylinder structures are formed, which finally relax to 
a vague bulge in 10® yr. (b) Relaxation is dominated 
by self-gravity of the shell under the infiuence of the 
disk potential. An expanding torus appears and a 
few clumps are formed in this torus near the disk 
plane. Finally, a compact high density core with a 
diffuse halo is formed by merging of these clumps. 

2. During the relaxation of the shell, the disk thick- 
ens. Non-axisymmetric structure (such as spirals) 
disappears due to the heating of disk stars by the 
gravitational interaction between the disk and shell. 
This heating is more effective in dynamically cold 
shell models, in which an oval density wave expands 
toward the outer region. In the relaxation process of 
the dynamically cold shell, a few clumps grow and 
heat up disk stars. 

3. Angular momentum of the disk is transferred to the 
shell stars during relaxation by dynamical friction 
between the clumps consisting of shell stars and disk 
stars. As a result, the final bulge begins to rotate. 
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although the shell initially had no rotation. Mean 
rotational velocity of the bulge is ~ 50 km s~^ at 
1 kpc, and the spread of the line-of-site velocity in 
their position-velocity map decreases with the ra- 
dius. These results are consistent with recent ob- 
servations of the bulge of our Galaxy (Nakada et al. 
1993; Minniti et al. 1992). 

5.2. The ongm of SO galaxies 

Formation of SO galaxies has been a controversial prob- 
lem. There appears to exist two scenarios for the forma- 
tion of SO galaxies. (A) The properties of SO, such as 
large bulge/disk ratio, were determined by the condition 
at galaxy formation. According to Larson (1976), the 
bulge/disk ratio depends mostly on a star formation rate 
during initial collapse of a proto-galaxy and how it varies 
with time: The spherical component is formed in an early 
phase of rapid formation of stars. In this mechanism, 
the large bulge/disk ratio of SO galaxies is caused by a 
high star formation rate in a high density proto-galactic 
core. (B) The large bulge/disk ratio was produced by 
the interaction between galaxies and their environment 
at the early stage of galaxy formation: This scenario in- 
cludes stripping of gaseous envelope of gas-rich galax- 
ies at collapse of a cluster of galaxies (Larson, Tinsley, 
Caldwell 1980), merging of proto-galaxies (Silk, Norman 
1981), accretion of small satellites (Pfenniger 1993) or 
starbursts induced by galaxy-galaxy interactions (Sofue, 
Habe 1992). 

We also suggest that the scenario of bulge formation 
following a starburst-induced stellar shell and its relax- 
ation is particularly important for producing a large 
bulge/disk ratio. Our simulations end up of looking like 
SO galaxies (e.g. de Carvalho, da Costa 1984; Hamabe, 
Wakamatsu 1989). Figures 11(a) and (b) are density con- 
tours of disk component projected onto the x-z plane for 
model S2 and S3 at t = 2.0, respectively. Figures 12(a) 
and (b) are the same, but for both disk and bulge com- 
ponents. These figures show that, while the disk compo- 
nents of the two models are quite different, the composite 
isophotes resemble those of SO galaxies. In the obser- 
vational sequence SOa - SOc (van den Bergh 1976) the 
bulge is most compact in SOc galaxies. Our models indi- 
cate that compact bulges are associated with thick disks, 
and vice versa. Namely, they suggest that the more mas- 
sive is the bulge compared to the disk, the thinner is the 
disk compared to the disk radius. In fact, we have found 
that the observed disk-to-bulge mass ratio appears to be 
inversely proportional to the scale thickness-to- scale ra- 
dius ratio of the disk (Wada, Habe, Sofue, Taniguchi, in 
preparation). 

Our model of SO formation is consistent with the key 
observational properties. There is no spiral in the disk 
because of interaction with the shell stars. A hot X- 



ray halo (Canizares, Fabbiano, Trinchieri 1987; Awaki 
et al.l991) originates in the superwind phase (Tomisaka, 
Ikeuchi 1988). If their structure is due to a process trig- 
gered by galaxy-galaxy interactions, the presence of SOs 
would increase with local number density of the galaxies, 
which is indeed observed (the morphology density rela- 
tion: Dressier 1980). A galaxy evolves into a gas-poor 
disk system, as the gas is used up in the starburst phase. 

We simply assumed an initially thin stellar disk with 
no spheroidal stellar component for the live disk model in 
section 4.2. This model is based upon Larson's model for 
forming disk galaxy with a small bulge/disk ratio (Lar- 
son 1976). This disk dominant galaxy can be formed 
from the collapse of a rotating proto-galactic gas cloud 
before star formation. Even if a small bulge already ex- 
isted before the stellar shell formed, our results would 
not change significantly. A shell with small velocity dis- 
persion would continue to be strongly affected by the 
disk potential away from the central region. Interaction 
between a shell and a disk with a massive pre-existed 
spheroidal component would be an interesting subject. 

The initial velocity dispersion of the shell stars is one 
of important parameters for dynamical evolution of the 
system. It would be generated by the Rayleigh- Taylor 
instability of the gas shell, and by the time difference of 
star formation in the decelerated gas shell, in addition to 
the star formation process itself (i.e. gravitational and 
the thermal instability of the gas shell). However, it is 
hard to determine the initial velocity dispersion of the 
shell stars either observationally or theoretically. In or- 
der to investigate this process theoretically, two or three 
dimensional hydrodynamical simulations of shell forma- 
tion by superwind taking into account self-gravity and 
star formation of gas would be necessary. 

There is a possibility that newly formed stars exist far 
above a disk of nearby starburst galaxies. For example, 
the central region of the typical starburst galaxy M82 is 
associated with a large-scale outfiow of dense molecular 
gas (Nakai et al 1987). It would be expected that this 
molecular cylinder is a site of star formation. In fact, ion- 
ized gaseous filaments are associated with dark filaments 
of M82 (Lynds and Sandage 1963), and UV radiation 
from newly born OB stars would be a possible excitation 
source. Detailed observations of the starburst regions 
will bring us important information about kinematics of 
the newly formed stars outside the disk. 

Some authors proposed that the bulge of some galaxies 
may be accumulated from the merger of satellites, such 
as globular clusters or dwarf irregular galaxies (Schweizer 
and Seitzer 1988). From a dynamical point of view, this 
model is basically the same as the collapse of the shell 
stars in our model. However, there is a serious prob- 
lem in this scenario. Stars in dwarf galaxies are metal 
poorer compared to more massive galaxies. This is be- 
cause the low star formation rate in dwarf galaxies due 
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to tlieir low ISM density, and mass of the dwarf galaxies 
is so small that they cannot confine metal enriched gases 
ejected from stars. On the other hand, bulge stars in 
our Galaxy are metal rich (Rich 1988). Chemical evolu- 
tion model of the bulge of our Galaxy suggests that the 
bulge was formed in the initial 10® yr; the IMF should 
have been more gentle than the Salpeter's IMF; and the 
time scale of star formation should be less than 10® yr 
(Arimoto and Yoshii 1987). This means that bulge stars 
were made from gases contaminated by gas originated 
in the primeval starbursts. Population synthesis simu- 
lations have also shown that a star formation history in 
a bulge is the same as that in the bulge of our Galaxy 
(Jablonka, Arimoto 1992). Consequently, recent merger 
events cannot produce the bulge of our Galaxy. On the 
other hand, in our starburst model, the gases in the ini- 
tial gas shell are mixed with metal rich superwind caused 
by the nuclear starburst, and it is natural that the bulge 
stars formed from the shell stars are already metal rich. 
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Figure captions 

Figure 1. Evolution for model C projected onto the 
x-z plane. Time is shown in the upper right-hand 
corner in unit of lO^yr. All frames are 20 kpc x 20 
kpc. 

Figure 2. Evolution for model B projected onto the x- 
y plane (left panels) and onto the y-z plane (right 
panels). All frame are 20 kpc x 20 kpc. 

Figure 3. Evolution for model A2 projected onto the 
x-y plane (left panels) and onto the x-z plane (right 
panels). All frames are 20 kpc x 20 kpc. 

Figure 4. M^/Md - a^^ dia gram, where Mg, M<j, and 
(Tg are masses of the shell and the disk, and the ini- 
tial velocity dispersion of the shell stars, respectively. 
The two horizontal lines, (T<j, derived from the z- 
component of the Jeans equation for a Miyamoto - 
Nagai potential disk at i? = 1.0 and at i? = 10.0. 
Unit of velocity is 10 km s~^. Relaxation types of 
models Ai , A2 , A3 , B, and C are classified into 
Type I, II, and III (see Section 4.1). 

Figure 5. Evolution of model S2 projected onto the 
x-y plane for its disk stars only {center), for shell 
stars only (left), and for both shell and disk stars on 
the x-z plane {right). Times are shown in the upper 
right-hand corner of the right panels. All frames are 
24 kpc X 24 kpc. 

Figure 6. Same as Fig. 5, but for model S3. 

Figure 7. Time variation of mean scale-height of the 
disks of models SI, S2, S3, NS. Unit of the scale- 
height is kpc, and that of time is 10® yr. Model NS 
is a comparison model as described in Section 4.2, 
in which the shell was added to the stellar disk at 
time = 0. The open circles represent a scale height 
of a disk without the shell. 

Figure 8. Evolutions of mean vertical velocity disper- 
sion of disk stars for models S2 and S3. 

Figure 9. Time variation of specific angular momentum 
normalized by the initial value around the z-axis for 
disk stars (left vertical axis) and for shell stars (right 
vertical axis) of models SI, S2 and S3. 

Figure 10. Density profiles of the final shell stars for 
models S2 and S3. Lines of p cx and p oc r~^ are 
also shown . Both bulges have a core (p cx r~^)-halo 
(p oc r~^) structure. The dynamically colder shell 
model(S2) has a smaller size and higher density core, 
comparing to S3 model. 



Figure 11. (a) Surface density contours of edge-on view 
only for the disk stars of model S2. (b) Same as (a) 
but for model S3. Size of the frames are 24 kpc. 

Figure 12. Same as Fig. 12, but for both disk and bulge 
stars. 



